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CN ■ Abstract 

On . 

We discuss the propagation of wave packets through interacting envi- 
Q_i| ronments. Such environments generally modify the dispersion relation or 

shape of the wave function. To study such effects in detail, we define the 
distribution function Px(T), which describes the arrival time T of a packet 
at a detector located at point X. We calculate PxiT) for wave packets 
traveling through a tunneling barrier and find that our results actually ex- 
plain recent experiments. We compare our results with Nelson's stochastic 
interpretation of quantum mechanics and resolve a paradox previously ap- 
%— i ' parent in Nelson's viewpoint about the tunneling time. 

1 Introduction 

We are interested in the behavior of quantum particles, that is, wave packets 
propagating through interacting environment. In general, there are two types 
of environment. One is the ordinary medium (plasma, dielectric, etc.) which 
consists of "matter" [1-4]. The other is the nontrivial structure of the vacuum 
due to field theoretical fluctuations p or effects of quantum gravity || In 
both cases, the presence of such environments will modify the dispersion relation 
of particles, E = f(p), or modify the shape of the wave packet. Observation of 
the arrival time of particles through such environments is a way to see the effects 
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of these modifications. Recently, these effects have been tested in two fields, 
astrophysics and quantum optics. The first is the observation of arrival times of 
photons from distant astrophysical sources such as 7-ray bursters. Several models 
of quantum gravity suggest that the velocity of light has an effective energy 
dependence due to the modified dispersion relation induced by the nontrivial 
structure of space-time at distances comparable to the Planck length. To confirm 
this effect, it is necessary to observe a certain difference of the arrival time of 
photons with different energies, and 7-ray bursters work for this purpose || . As 
a result, a lower bound on the energy scale of quantum gravity is obtained ||. The 
second recent test is observation of tunneling of photons. Chiao and co-workers 
constructed an elaborate stadium for the race between photons propagating in 
the vacuum and through an optical barrier, and measured their arrival times 
[|l], 0. They found that the photon tunneling through the barrier arrived at the 
goal earlier than the other photon traveling in the vacuum. Although this result 
implies superluminal velocity of the tunneling photon, it does not mean causality 
violation, because in this case the group velocity itself does not transport any 
information at all. The apparent superluminality results from reshaping of wave 
packets while tunneling. Similar phenomena can be found in absorbing media ||. 
Anyway, in both experiments, measurement of the arrival time of wave packets 
plays an essential role. 

However, there is no clear definition of arrival time in quantum mechanics. 
This has its root in the well-known fact that time is not an operator but a pa- 
rameter in quantum mechanics. Though many authors have attempted to define 
an operator of arrival time and construct its eigenstates, a satisfactory formula- 
tion has not yet been obtained [10-25]. In this article we define a distribution 
function Px{T), which describes the arrival time of packets at a detector located 
at point X. In terms of Px(T), we can compute a mean arrival time (T) x . Of 
course we assume an ideal detector and our definition of Px(T) might not exactly 
correspond to the physical measurement process. However, concrete calculation 
of Px{T) shows us clearly the dynamical properties of propagation of packets 
through interacting environments. 

We investigate the arrival time distribution Px{T) numerically for nonrela- 
tivistic massive particles traveling through a potential barrier in one space di- 
mension, that is, tunneling processes. This might be a simple model for the 
experiment of Chiao and co-workers. In this case the existence of a potential 
barrier V(x) causes reflection and transmission of packets; therefore the behavior 
of Px(T) will be highly nontrivial, depending on various parameters. How to deal 
with time in tunneling processes is also known as the tunneling time problem. 
The problem arises from the paradox that a particle under a potential greater 
than the particle's energy seems to move with a purely imaginary velocity. In 
recent developments of nanotechnology, the study of the tunneling time has great 
significance because it might enable us to estimate the response time of nanode- 
Z7|| . Various approaches to the tunneling time have been proposed by many 
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authors [27-36] ; however, it seems difficult to define it uniquely Q Therefore we 
need to define effective tunneling times for each system and each purpose. We 
have no intention of wrestling with the general theory of tunneling time now; 
therefore, we restrict ourselves to analyzing the time of appearance of the packet 
in the exit of the potential barrier and how it moves after that. These two notions 
determining the arrival time difference have usually been confused. In this article 
we will distinguish them clearly. 

Finally we consider the real-time stochastic interpretation of quantum me- 
chanics introduced by Nelson [p7|| . Since it utilizes the real-time trajectories of 
quantum particles as sample paths, we can construct an appropriate time distri- 
bution from ensemble of sample paths. This is why Nelson's approach is expected 
to be effective for time problems in quantum mechanics. In particular, it is in- 
teresting to attack the tunneling time problem from this approach because we 
can trace the particle's real-time motion even under the tunneling potential. Ac- 
tually it has been found that the tunneling particle "hesitates" in front of the 
barrier [|38|] . This property seems paradoxical because it implies that the particle 
tunneling through the barrier should always be delayed compared with the free 
one due to this hesitation and it seems contradictory to the advancement of the 
peak of the wave packet as seen in the experiment of Chiao and co-workers. Is it 
a real paradox? 

It is clear that Nelson's approach can reproduce any physical quantities of the 
usual quantum mechanics by averaging them about the sample path ensemble. 
However, there is no reason that any "observables" classically defined in Nel- 
son's stochastic procedures should have corresponding quantities in the standard 
quantum mechanics. We will compute the arrival time distribution in Nelson's 
approach and compare it with our Px{T). Then we clarify the real physical mean- 
ing of the "hesitation" and show that there is no paradox at all. Furthermore, 
we mention that Nelson's interpretation can explain the characteristic behavior 
of (T) x for tunneling particles very well. 

2 Definition of the arrival time distribution 

First we will briefly review previous attempts to define a time of arrival op- 
erator and their difficulties. In the 1960s, Aharonov and Bohm quantized the 
representation of the classical arrival time for the free particle at a point X = 



x » m/,1 1,\ 

T = —m— — > T = x — I — x . (1) 

p 2 \ p p J 



1 "The systematic projector approach" has been proposed as a unifying theory of the various 
times proposed so far [B3|. 
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Here x and p are the initial position and momentum, respectively, where we 
work in the Heisenberg picture. Because T satisfies T,H = i%, it seems a 

good definition. We construct its eigenstates T\T) —T\T). Q However, these 
eigenstates turn out to be not orthogonal, 

(p\T) oc [6{p) + i6{-p)]^ 2T l 2m \ (2) 
(T|T') cc *(T - TO - £P— j^, (3) 

where P represents Cauchy's principal value. That is, T is not Hermitian. The 
origin of difficulty is the singular behavior of T at p = 0. Recently the regulariza- 
tion of T with an infrared momentum cut off [15| and an interpretation by means 



of the positive-operator- valued measure were proposed However, the valid- 
ity of this procedure is not clear ]17| [18|] . In the first place, there is no one-to-one 
correspondence between the operator representation in quantum theory and the 
classical representation, and it becomes more complicated for interacting cases 
[19-22]. 

Now we will not insist on defining an arrival time operator; rather, we try to 
construct an arrival time distribution directly. We suppose that there is a detector 
on the path along the motion of wave packets and it counts the particle according 
to the value of the wave function ip(X,t) at every time t = T. Supposing the 
detector is ideal, we directly define the arrival time distribution Px{T) from 
^(X,T), 

P x {T)dT= r L x{T)dT , p x (T) = \ij(X,T)\ 2 . (4) 



roc 

/ dT Px (T) 
Jo 



Although Eq. (f|) looks like a trivial definition in our picture, we will derive 
it, clarifying our system setup and assumptions. We consider a system consisting 
of a particle and a detector located at x — X. If there is no interaction between 
them, the system Hamiltonian Hq and the system state |\&) are given by 

H = H p ® 1 + 1 ® H D , (5) 
I*) = (6) 

where H p is the particle Hamiltonian, is the particle state, and similarly H^ 
and \D) are those of the detector. We define the total Hamiltonian H by adding 
the interaction Hamiltonian H\ between the particle and the detector, 

H = H + H h Hj = gV p (x) ® V B . (7) 



2 In order to obtain a complete set, one needs two eigenstates \T, ±) for every value of T lli 
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For simplicity, we consider a detector whose state consists essentially of two 
components, 



II) 




1 



1 




, reacted. 



, unreacted, |f) 

Corresponding to this representation, we set the interaction potentials 

V p (x) = 5(x-X), 

V D = 



(8) 



1 

1 



IT)- 



(9) 
(10) 

This choice of Vd should be 



which induces a transition ||) =>- Vb |!) 
meaningful only in the first order of g. 

Now we consider the time evolution of the system from t — to T. We 
prepare the initial state \D(0)) = \l) and evaluate a quantity Rx(T), which is 
the probability that the state \D(t)) is found to be | j) when t = T. From Rx(T), 
we get Px(T)AT, which is the probability that the transition \l) =>- ||) occurs 
in a time interval [T, T + AT], that is, 

Px(T)AT = R X {T + AT) - R X {T) P x {T)dT = dR x (T) . (11) 

Next we evaluate Rx(T) in terms of the particle wave function ip(x,t). We 
now assume that the detector reacts only once incoherently, and therefore we 
calculate only in the first order of g. Adopting the interaction picture, the time 
evolution of the state can be represented as follows in the first order of g: 

= Te- l/n /o Td '^Pi(^)®vbiW |^(0)) l( g) |T>(0))j 
~ 1^(0)): ®|T>(0))i 



T 



~h Jo dt9V ^ x > ^ IV;(0)) i Vm ^ |jD(0)) i 



(12) 



where T represents the time ordered product. |\&(0))j is the undetected state and 
|^(T))j is the detected state, which is written in the Schrodinger picture as 



|¥(T)> = 



T 



dt 



e -iH p (T-t)/H gVp ^ e -iH p t/H^ |^ Q ^ 

= [ T dt \^(T;t))®\D(T;t)), 

n Jo 

where we introduced 

\1>(T;t)) = 



iH p (T-t)/n 



gV p (x)e 



-iH p t/ti 



\D(T;t)) = \e- iH ^ T ' t)/h V D e- tHot/h ]\D(Q)) 



(13) 

(14) 
(15) 



5 



We obtain Rx{T) in terms of the norm of the detected state, 



R X {T) = (*(T)\*{T)), (16) 

under our approximation of weak coupling. Now we apply a macroscopic deco- 
herence condition, 

(D(T;t 1 )\D(T',t 2 )) = 6(t 1 -h). (17) 

This means that the states reacted at different times are orthogonal to each other, 
that is, once the detection process occurs, the total state effectively loses its 
coherence and looks like a mixed state. Of course it is not possible to satisfy this 
condition by working in the two-dimensional Hilbert space in Eq. (P) . We should 
describe the detector by means of an infinite-dimensional Hilbert space to realize 
decoherence effectively ^6|. However, we can avoid this assumption by "switching 
on" the interaction Hamiltonian during different small time intervals in repeated 
experiments, instead of using the finite time interval [0, T] and differentiating 
with respect to T. 

Under this condition, the evaluation of Rx{T) and Px{T) is straightforward 
as follows: 

Rx(T) =Wrj\W)) = i f T dh [ T dt 2 (^(T;t 1 )\i;(T;t 2 ))8(t 1 -t2) 

n Jo Jo 

= 1 [ T dt ^(T;t)\i>(T;t)) , (18) 
n Jo 



d „ _ d 



Px(T) = —R x (T) = — (*(r)|*(T)) = ^^(T;T)|^(r;T)> 

= (m \e lH ^ h V^x)V p (x)e-* H ^ h \ ^(0) 
= ^(0)|^(X,T)| 2 , (19) 

where at the last step we inserted the complete set (/ dx\x) (x\) three times. 
Although the divergent 5(0) seems to break the validity of our formulation, we 
can remove this singularity by replacing the 5 function in Eq. @ with a smeared 
function. We normalize the right hand side of Eq. fll9| ) to get our expression 
for the arrival time distribution Px(T). Using Px{T) we define the mean arrival 
time (T) x , 

rod 

(T) x = / TP x (T)dT. (20) 

•/ 

Because Px{T) and (T) x have simple and general expressions, we can calculate 
them easily even for interacting cases. Our Px{T) is often called the "presence 
time distribution" because of its behavior in the classical limit |2(| . In order to 
avoid confusion, we should make clear that the distribution Px{T) may not be 
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interpreted as the probability distribution of a quantum mechanical time observ- 
able. It is an effective distribution describing a "relative probability." 

Of course, our definition of arrival time distribution Eq. (|j) is not a unique 
one. Considering a different system setup, some people have proposed a definition 
using the current Jx{T) instead of px{T) = |^| 2 [ P"l|j , 

P X (T)dT= /j^ dT , M T) = hm( r d /] . (21) 
/ dTJ x (T) m V 9xJ x=x 

Jo 

This definition has the serious problem that Jx (T) can be negative in some cases, 
for example, detection before the potential barrier. Therefore we cannot identify 
P£(T) as a probability distribution. As for detection beyond the potential barrier 
as we discuss below, Jx{T) might effectively maintain positivity and actually the 
behavior of P X (T) is found to be similar to ours. 



3 Calculation of Px(T) for tunneling particles 

Now let us calculate Px(T) and (T) x for nonrelativistic massive particles trav- 
eling through a potential barrier V(x) in one dimension. This is a simple model 
of tunneling processes such as the experiment of Chiao and co-workers . Solving 
the time dependent Schrodinger equation with some initial conditions, we can 
get ip(x,t). Except for the free case it is difficult to solve the partial differential 
equation analytically, and therefore we solve it numerically. We now employ a 
discretization scheme known as the Crank-Nicholson method, which conserves 
the norm of ip(x,t) even with a finite discrete time step j|2]. We work with the 
units m = h = 1 and for the initial condition we prepare a Gaussian wave packet, 

ij)(x,0)=[ — ) e -(*-*o)V2^ e «D(*-*o) , 22 ) 
\na 2 J 

whose mean energy is (E) = + 1/4<t 2 , and we set a time independent 

square potential barrier V(x) on a section [0,d\. For simplicity, in this article 
we work with a unique initial packet. We fix the central wave number k = 2 
and in this unit we set the width of the initial packet in the configuration space 
a = 10(2//co) and the center of initial packet xo = —50(2/k ). All quantities that 
have a dimension of time are measured by the (4//cq) unit. Hereafter we will 
omit the units of the numerical values. We change two parameters of the barrier 
potential V(x): the width d and the height h, and also the detector location X. 

Let us begin by watching the motion of wave packets with an h = 1.1(E), 
d = 1.5 potential. The snapshots of the motion are shown in Fig. [l], in which 
the free packet motion is also shown for comparison. Both packets spread due to 
the dispersive properties that come from their own masses. The packet moving 
through the potential barrier experiences reflection and transmission and the 
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h=1 .1 <E> , d=1 .5 [2/ko] h=2<E> , d=4 [2/m 




Figure 1: Snapshots of the wave func- Figure 2: The arrival time distribution 
tion squared at various times for detection at X = 50[2/A;o] 

peak of the transmitted part will often advance compared to the free packet . It is 
usually explained that this is because the higher momentum components of the 
packet preferably go through the barrier and they propagate faster than the lower 
momentum parts due to their dispersive properties. That is, the advancement 
results from reshaping of the transmitted packet. However, tracing the peak of 
the packet is often difficult because near the barrier the peak cannot be clearly 
identified. Therefore we must use more well-defined quantities, Px{T) and (T) x . 

Analysis 1: Detection at X = 50 

Now let us calculate Px(T) and (T) x at X = 50 with an h = 2(E), d = 4 
potential. In Fig. ^|, the arrival time distributions Px{T) are plotted for the free 
and tunneling particles and the mean arrival times (T) x are shown by dashed 
lines. The remarkable feature of Px{T) is the stretched tail and the shift of the 
peak caused by spread of the packet. For the free case, (T) x = 50.13 is later than 
T = 50, which is expected from the group velocity of the free packet, and "the 
peak of Px{T) n = 49.94 is earlier than T = 50. It is also seen that, because of the 
packet's reshaping, Px{T) for the tunneling particle has a narrower shape than 
the free one and U (T) X for the tunneling particle" = 47.65 is earlier than the free 
one. However, it should be noted that only one detection far from the barrier 
cannot describe the dynamics of packets since we should discriminate effects in 
and out of the potential barrier. Therefore we investigate detection at various 
points. 

Analysis 2: Detection at various points 

When we try to give a definite answer to the so-called tunneling time problem, 
we might have to calculate the difference A = (T) d — (T) , since this problem 
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demands that we answer the question "How long does it take for the particle 
to tunnel across the barrier?" However the difference A does not make much 
sense because, as we can see in Fig. |], the shape of the packet is oscillating 
frequently at the entrance of the barrier and it is difficult to distinguish between 
the tunneling packet and the reflected one, that is, (T) is not a good physical 
quantity. On the other hand, the packet has a relatively clear shape at the exit 
of the barrier. Therefore we can analyze what time the packet will appear at the 
exit of the barrier and how it moves after that. 




We calculate (T) x for the exit of the barrier and several points after that: 
(T) x at X = d, 10,20,30,40,50, with d = 4, h = 0.5(E), 1.1(E), 2(E) barrier 
potentials (Fig. [|). We can see two remarkable features in this figure. The first 
is that for high barriers (T) d is earlier than in the free case, but it is later for 
low barriers. That is, it seems that the transmitted packet arrives at the barrier 
exit earlier than the free one for tunneling dominated cases. These are regarded 
as effects in the barrier. The second feature is that, after passing the barrier, 
the tunneling packet moves with a constant mean velocity larger than that of 
the corresponding free packet. This is an effect outside the barrier. These two 
types of effect are combined to cause nontrivial behaviors of the arrival time. For 
example, in the case of d = 4, h — 0.5(E), the tunneling packet arrives at the 
barrier exit X = 4 later than the free packet; however, after exiting the barrier, 
the tunneling packet catches up with the free one and overtakes it at X ~ 15. 
After all, it depends on X which arrives at X earlier, the tunneling or the free 
packet. 

We can see the second effect clearly in the Fourier transformed form of the 
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transmitted packet |33 

= (47RX 2 ) 1 / 4 1 ^ e - 2 (fc-fco) 2 /2| Tfe | e ^ e ,[fc(.-.o)-^ ; (23) 

where Tj, is the transmission amplitude and 9 is the phase. Using an analytically 
obtained \Tk\, the mean momentum k m can be calculated for the transmitted 
packet. Results for several potential conditions are shown in Fig. For "high" 
barriers, a wider barrier gives a larger mean momentum in the region d : [0,4]. 
This is because as d grows the \Tk\ support shifts to the higher momentum side. 
Therefore the statement "Higher momentum components of the packet preferably 
go through the barrier" applies indeed. This kind of "acceleration" effect is found 
in other areas of physics [ p8fl . 

Analysis 3: Detection at the barrier exit X = d 

To see the in-barrier effects more definitely, we calculate the difference between 
mean arrival times for the tunneling packet and the free one at the barrier exit 
X = d, 

AT = (T)* unncl - (T)J ee . (24) 

Results for the same potential conditions as in Fig. ^ are shown in Fig. [5[ At 
first we see that in the small d region AT is positive, that is, the tunneling packet 
gets behind the free one, for any potential height. However, as d increases, AT 
shows different behaviors according to the potential height. 

Roughly speaking, for a "low" barrier AT almost stays positive but for a 
"high" barrier AT becomes negative. The "high" barrier means that the tun- 
neling modes dominate in the transmitted packet. In the large d region, AT is 
negative, that is, the tunneling packet goes ahead of the free one for the tunneling 
dominated case. We also see a strange behavior where AT changes sign twice and 
finally becomes positive. The typical case in Fig. is the h = 1.1(E) barrier. 
We understand this effect as follows. For a very wide barrier, over-the-barrier 
modes dominate in thetransmitted packet (ou m = k^/2 > h). That is, as in the 
"low" barrier case, AT becomes positive again. 

In Fig. |5] we also plotted an analogous quantity AT^, calculated by the sta- 
tionary phase method. We define AT^, as follows: 

T( f = ~J~. > ( 25 ) 




AT V = — — — (d - x ) + r v , (26) 



where 6 is the phase shift of the transmitted wave defined in Eq. (^3|) and v g (-) 
are the group velocities v g (k m ) = (du/dk)\ k=k = k m , v g (k ) = k . In the 
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Figure 5: The difference between mean arrival times at the barrier exit X = d for 
the tunneling packet and the free one: AT = ^^ unnel _ (T)^ ree . AT V (shown by 
dots) are the same quantity calculated by the stationary phase method for each 
potential. 



ordinary tunneling time problem context is called the phase time. As seen in 
Fig. [|, although AT V has good agreement with our AT in the small d region, as 
d increases, the difference becomes clear for h (E) barriers. This is because for 
such barriers the momentum distribution e~°" ( fc ~ fc °) / 2 |Tfc| is no longer symmetric 
with respect to k m , and the packet's peak given by the stationary phase method 
loses physical significance. 

We would like to close this section by referring to the relationship between 
our results and the experiment by Chiao and co-workers, that is, tunneling of 
the massless photon. Of course our model does not describe the propagation 
of photons, and we now mention only the qualitative behavior. Because the 
energy of the photon in the vacuum is exactly proportional to its momentum, 
the group velocity of the photon after tunneling is constant c. Therefore we 
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get an X independent constant value of the difference AT — (T) x — (T) x 
at any X > d. Since their experimental setup is the tunneling dominated one, 
it may correspond to our model with high and medium wide barriers. Then 
our results are consistent with their experimental observation that the tunneling 
photon arrives earlier than the free photon. 



4 Nelson's stochastic interpretation 

Now we consider the stochastic interpretation of quantum mechanics intro- 
duced by Nelson. This approach interprets the motion of particles in quantum 
mechanics as "real-time" stochastic processes P?| . Nelson substituted the coor- 
dinate variable x(t) for a stochastic variable performing the Brownian motion in 
a certain drift force field. The time evolution of x(t) is described by the Ito-type 
stochastic differential equation, 

dx{t) =b(x(t),t)dt + dw(t), (27) 

where b(x, t) is the so-called drift term, given by the ordinary Schrodinger wave 
function ip(x, t) as 

h d 

b(x,t) = — — (Im + Re)lm/>(x,t). (28) 
m ox 

The Gaussian noise dw characterizes the stochastic behavior and should have the 
following statistical properties: 

h 

(dw(t)) = 0, (dw(t)dw(t)) = —dt. (29) 

m 

Starting with an initial distribution of x(0) we solve Eq. ( p7[ ) and obtain sample 
paths. Averaging a physical variable with these sample paths, we can calculate 
the expectation value for the ordinary probability distribution \ip(x, t)\ 2 . In this 
approach, we are able to observe "trajectories" of real-time motion of a parti- 
cle, that is, to describe the quantum mechanical time evolution by a classical 
stochastic process. 

Thus in Nelson's approach it may be possible to understand an imaginary- 
time process such as tunneling in real-time language. It was pointed out that 
the tunneling particle "hesitates" in front of the barrier as seen in Fig. [| |p8|| . 
This fact was understood to imply that the particle tunneling through the barrier 
should always be delayed compared with the free one because of this hesitation. 
Is it contradictory to our results? Nelson's approach can reproduce physical 
quantities in standard quantum mechanics, and there cannot be any conflict. 

Now we analyze the mean arrival time in Nelson's stochastic interpretation. 
One intuitive idea of defining the arrival time for a sample path is to measure 
the time for a path to reach a detecting point for the first time: "the first time 
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h=2<E> , d=1 [2/M 




counting scheme" 39| . However, this notion has no counterpart in the physical 
quantities of standard quantum mechanics. We have to work with the probability 
of existence of paths at a point (or a section) at some definite time. The difference 
between these two notions is that the latter counts the possibility of a path going 
beyond the point and coming back to it at the measuring time. 
We define a probability function p x (T), 

P mdx=^i, (so 

where N is the total number of sample paths and n(X, T) is the number of sample 
paths that exist in [X, X + dx] at time T. As stressed before, we will count the 
number of paths passing a target point over and over again, i.e., we now employ 
"the multiple counting scheme." With this scheme, we define the arrival time 
distribution P X (T) and the mean arrival time (T) x of the particle in Nelson's 
stochastic interpretation, 

(T) dT = r P J {T)dT , (31) 
dTpx (T) 



o 



roc 

(T) x = / TP^{T)dT. (32) 



o 



We calculate P X (T) and (T) x with an h = 2 (E), d = 1 barrier by solving Eq. 
to get iV = 10 6 sample paths. The result is shown in Fig. [7[ The distribution 
Px(T) agrees with Px{T) very well; therefore {T) x agrees with {T) x . The 
distribution given by Nelson's approach exactly reproduces our previous results, 
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h=2<E> , d=0.5 [2/M 



h=2<E> , d=1 [2/M 
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Figure 8: Sample paths for AT > Figure 9: Sample paths for AT < 



just as expected. Of course, if we employ the first time counting scheme, -Px(T) 
shifts to an earlier time region and therefore (T)^- is smaller than (T) x . 

Then, we should answer the paradoxical question, "Why does a hesitating 
particle arrive earlier than the free one?" To answer this question, let us compare 
the two cases of AT > and AT < 0. First we show typical sample paths for 
two cases, AT > (AT = 0.084), h = 2 (E), d = 0.5 in Fig. | and AT < 
(AT = -0.138), h = 2 (E), d = 1 in Fig. g. In both figures we also plot the 
average position of the free sample paths (x) free . We should pay attention to the 
point (x,t) = (0,25) because (x) lree arrives in x = at t = 25. The two figures, 
Fig. § and Fig. |9|, make a remarkable contrast. That is, although in Fig. |8| even 
the paths that arrive at x = later than t = 25 can pass through the barrier, 
in Fig. |9| essentially only the paths that arrived at x = earlier than t = 25 
can go through it. The "hesitation" property is seen in both cases. In Fig. g, 
however, even with "hesitation," the averaged tunneling path can appear at the 
barrier exit x = d earlier than the averaged free path because the tunneling path 
arrived at x = much earlier than (x) {vee . This is the key to the mystery between 
hesitation and advancement. 

Well, why do the tunneling paths conduct themselves in such a strange way? 
In the first place, why does hesitation occur? The reason is hidden in the time 
dependence of the drift velocity b(x,t). We show b(x,t) for the same conditions 
discussed above, especially near the potential barrier (Fig. [10] and Fig. |il"|) . 



In the foreground of the barrier, according to the interference of the incident 
packet and the reflected packet, b(x,t) oscillates frequently and becomes null 
many times. Especially near the barrier entrance x = 0, b(x, t) changes from 
positive to negative, where the particle is "trapped." These effects cause the 
path's hesitation. 

At earlier times, b(x, t) is almost always positive value but at later times, 
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Figure 10: Drift velocity for AT > Figure 11: Drift velocity for AT < 



it becomes almost always negative. In Fig. [11], this tendency is extreme and 
realization of the tunneling path is much rarer than in Fig. [H]. This is the reason 
why the early arrived paths tend to pass the barrier more easily. After all, there 
is no inconsistency between our results and the hesitation behavior in Nelson's 
interpretation. 

Furthermore, Nelson's interpretation provides us an intuitive explanation of 
our results. Let us consider the high potential barrier case. It is important that 
every transmitted path hesitates to some extent. In the small d region, because 
of the high transmission rate, even a path arriving at x ~ relatively late can 
pass the barrier, and as a result we find AT > 0. As d increases, the transmission 
rate becomes lower and only the paths arriving at x ~ earlier can penetrate 
the barrier, and as a result we find AT < 0. Finally, as d becomes very large, the 
paths arriving at x ~ very early hesitate there for a very long time; therefore 
AT becomes positive again. 

Of course, we must recall that the "path" in Nelson's view never corresponds 
to a real particle in ordinary quantum mechanics, and the explanation we gave 
above is just an interpretation. The same is true for an interpretation by the 
Bohm trajectory [|l!| It may be interesting to regard the path as a physical 
one and to calculate various quantities that cannot be calculated in ordinary 
quantum mechanics (the tunneling time A N = (T) d — (T)q , quantities calculated 
in the first time counting scheme, etc.). Although these attempts may give us 
deeper insights into quantum dynamics, the validity and significance of them have 
not been argued much so far [||| fill . 
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5 Summary 
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Supposing an ideal detector, we defined simple expressions for the arrival time 
distribution P x and the mean arrival time (T) x , and applied them to analysis of 
wave packet tunneling. We defined AT = (7^* unnel _ (T)^ ree and calculated it for 
various barrier conditions and showed the barrier effects clearly. In the small d 
region, AT is always positive, but as d increases, AT > for the over-the-barrier 
case and AT < for the tunneling case. After tunneling, the packet usually 
moves faster than the free one because it preferentially consists of the higher 
momentum modes of the incident packet. The barrier works as an acceleration 
filter in a sense. We also clarified that the stationary phase method gives a good 
approximation to our results, particularly in the small d region. 

We also confirmed that the stochastic interpretation introduced by Nelson 
reproduces our results. Furthermore, we clarified how the "hesitation" of the 
tunneling paths in Nelson's picture is consistent with the advancement of the 
tunneling packet. The key observation is that the paths arriving at the barrier 
earlier than the free mean paths tend to penetrate the barrier more easily. We 
pointed out that this property can be explained by the time dependence of the 
drift velocity b(x, t) and found that the behavior of AT is intuitively understand- 
able with Nelson's language. 
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